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Future interferometric gravitational wave detectors will make use of the coupling between shot 
noise and radiation pressure noise that produces a squeezed output for the quantum noise at the dark- 
port of the interferometer allowing these interferometers to operate with a sensitivity that exceeds 
the standard quantum limit at certain frequencies. The ability of the detector to take advantage of 
this squeezed output state depends in part on the readout system used. With the phase modulated 
laser input planned for these interferometers the quantum-noise limited sensitivity is slightly better 
when homodyne (rather than heterodyne) detection is used. We show that by modifying the laser 
input spectrum for these interferometers, readout with heterodyne detection can be made to have 
\Q • the same quantum-noise limited sensitivity as readout with homodyne detection. 
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The prospect of operating interferometric gravitational wave detectors for quantum non-demolition measurements 
| has lead to investigation into the details of the interferometer read-out scheme 0, • Two readout schemes relevant 
■ for gravitational wave detection are heterodyne readout (often called RF readout) and homodyne readout (often 
called DC readout). The readout scheme converts the field at the output port of the interferometer into a signal. 
In heterodyne readout a local oscillator field or fields that have a frequency difference from the field being measured 
interfere with the output field to produce an intensity that has a component which beats at the frequency difference 
of the fields. The phase of the beat note is a measure of the phase of the output field. In homodyne detection a 
local oscillator field at the same frequency as the field being measured interferes with the output field to produce 
, constructive or destructive interference depending on the relative phase of the output field. 

A principle issue for interferometric gravitational wave detectors is the difference in the quantum-noise-limited 
sensitivity of an interferometer with homodyne readout versus heterodyne readout. Buonanno, Chen and Mavalvala^ 
show that an interferometer with homodyne readout has better quantum-noise limited sensitivity than with heterodyne 
^Jy readout. Although other technical tradeoffs exist between homodyne and heterodyne detection we show that with 
a change to the laser input field assumed by Buonanno, Chen and Mavalvala it is possible for an interferometer 
with heterodyne detection to reach the same quantum-noise-limited sensitivity that can be achieved with homodyne 
detection. 
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II. INPUT SPECTRUM 



The conventional laser input spectrum for interferometric gravitational wave detectors is shown in the left-hand 
panel of figure ^ It consists of a high-power carrier at frequency loq which is resonant in the arms and interferes 
destructively at the dark-port of the interferometer, along with a pair of sidebands at ujo ± w TO that do not resonate 
in the arms and are made to transmit efficiently to the dark-port of the interferometer via an asymmetry in the arm 
lengths. Depending on the implementation there may be additional pairs of sidebands at other frequencies that are 
used for alignment sensing and control, which are not relevant for the discussion of the readout noise and which we 
shall ignore in our analysis. In heterodyne detection one or both of the sidebands at loq ± uj m are transmitted to the 
dark port of the interferometer where they interfere with the gravitational wave signal sidebands of the carrier field at 
frequency u>q ± Q gw , where Q gw <C u> m is the gravitational wave frequency which is small compared to the modulation 
frequency. The sidebands serve as the local oscillator for heterodyne detection of signals in a frequency band centered 
on the carrier frequency. With this input laser spectrum the signal sidebands at frequencies near ujq are detected, 
along with the vacuum fluctuations (quantum noise) that have frequencies in a band centered on cjq, as well as the 
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FIG. 1: Left: The input spectrum commonly used in gravitational wave interferometers. It has a high power component 
at frequency ujq which resonates in the arms of the interferometer and low power modulation sidebands at u>o ± <^m which 
transmit to the dark port and are used as a local oscillator for heterodyne detection. We refer to this as the "conventional 
input spectrum". Right: A "suppressed carrier input spectrum" with a low power component at frequency ujq which transmit 
to the dark port to be used as a local oscillator for heterodyne detection, and high power modulation sidebands at luq ± uo m 
which resonate in the arms of the interferometer 

fluctuations that have a frequency in either of the frequency bands centered on loq ± 2u> m . Since the spectral region 
around loq ± 2uj m contributes noise when it mixes with the local oscillator fields at loq ± <-o m and is demodulated, but 
does not contribute any signal, it degrades the quantum-noise-limited sensitivity of the interferometer. By contrast 
with homodyne detection a field at the carrier frequency luq is used as the local oscillator so that only the signal and 
quantum noise in the spectral region around uj$ contribute. Thus with homodyne detection there is less quantum-noise 
added to the signal. 

The carrier-suppressed input spectrum shown in the right-hand panel of Figure ^ leads to a difference in the way 
vacuum fluctuations affect the sensitivity. It consists of a low powered local oscillator at luq and two high-powered 
sidebands at ujq ± oj m . The sidebands are assumed to be resonant in the arm cavities while the local oscillator at 
frequency loq may or may not be resonant in the arm cavities depending on the implementation. Like with the 
conventional input spectrum, the Michclson interferometer's interference condition is set so that the high power 
components that resonate in the arm cavities interfere destructively at the dark-port while the low-power component 
that is used as a local oscillator is transmitted to the dark port of the interferometer. For comparison to the case of 
heterodyne readout with the conventional input spectrum or homodyne readout, we shall assume the same total power 
in the input spectrum is the same, so that the sidebands of our input spectrum each have half the power of the carrier 
in the conventional input spectrum (and the low power local oscillator components are assumed to have negligible 
power). With the suppressed-carrier input spectrum the local oscillator detects the gravitational wave signal sidebands 
in the frequency bands centered on loq ± uj m as well as the vacuum fluctuations at these same frequencies. Unlike 
heterodyne detection with the conventional input spectrum, there is no additional quantum noise from frequency 
regions that do not contribute to the readout of the gravitational wave signal. To analyze the quantum-noise limited 
sensitivity of a configuration with this supressed-carrier input spectrum we must consider the input-output relations 
for quantum noise in frequency bands centered on both the upper and lower laser sideband frequencies. We extend 
the analysis of the input-output relations for a single frequency range done by Buononno and Chen Q in which the 
quantum noise is represented by a 2-element vector and the interferometer is represented by a 2x2 matrix by describing 
the quantum noise with a 4-element vector and representing the interferometer by a 4x4 matrix. 



We consider vacuum fluctuations at the dark port of the interferometer as the source for the quantum noise that 
is read out. Specifically we look at the components of the vacuum fluctuations at frequencies loq ± oj m which are 
in-phase and quadrature-phase with the high-power sideband at each of those frequencies. We call the amplitude of 
the vacuum fluctuations that have a frequency in a band centered on the lower laser sideband frequency loq — uj m 
and are are in-phase and quadrature phase with this sideband a\ and 0,2 respectively. We let 0,3 and 0,4 represent the 
components of the vacuum fluctuations at frequencies in a band centered on the upper sideband frequency of luq + uj m 
that are in-phase and quadrature phase respectively relative to this sideband. With this notation, we consider the 
fields that contribute to the signal readout by heterodyne detection as a 4-element vector. The interferoemter can 



III. INPUT-OUTPUT RELATIONS 
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then be described as a 4 x 4 matrix that transforms the (unsqueezed) input quantum noise into a squeezed state 
at the output. By considering the action of a gravitational wave induced motion of the arm cavity mirrors on the 
output fields we relate the quantum-noise after demodulation to an equivalent signal magnitude to determine the 
quantum-noise limited sensitivity. 



^ arm cavity J 



signal extraction mirror 



FIG. 2: Diagram of a resonant sideband extraction interferometer. The operating point of the interferometer is set such that 
the laser power is reflected back towards the laser, and vacuum fluctuations from the dark port a; are reflected to the dark 
port bi 

We start with the input output relations for a conventional interferometer calculated by Kimble et al (eqaution 16 
ofRef. 0) 

di = Adi = cie 2t/3 

d 2 = Ad 2 + y/^JLe^, Ad 2 = (c 2 - KCi)e 20 (1) 

n a ql 

with the notation change of b — > d and a —> c to represent the internal fields of the interferometer configuration shown 
in figure [3 We extend these equations to account for the additional frequency components of the light circulating in 
the interfeormeter arms as follows. First these input-output relations which describe the fields at one frequency of 
interest (i.e. that of the lower frequency sideband, ujq — uj m ) are used to describe a second frequency of interest (i.e. 
that of the upper frequency sideband, u>q +u) m ) by making the replacements d\ — > d%, c\ — ► C3, d 2 —>■ d^ and c 2 — > C4. 
Next we add a term which accounts for a cross-coupling between these frequencies bands as follows: Differential motion 
of the arm cavity mirrors couples fields out of the interferometer into the dark port with an efficiency proportional to 
the magnitude of the field in the arms and represented by i/k. Analogous to equation [18] of Ref. [lj we set 

_ (I ± /I sq i)2^ 

* ft 2 ( 7 2 + ft 2 ) () 

Where I± is the intensity of the upper and lower sideband, I sq i — m 4t ^ 7 is the intensity needed to reach the standard 
quantum limit at the central frequency ujq for an interferometer with arm lengths L and mirrors suspended as free 
masses with mass m. f2 is the (gravitational wave induced) sideband frequency measured relative to the laser sideband 
at ujo ± Mm- 27 = j£ is the arm cavities bandwidth when the mirror that couples light into and out of the cavity has 
a power reflectivity of T. The in-phase vacuum fluctuations entering the dark port interfere with the laser power in 
the arms to produce radiation pressure noise that drives the mirrors differentially by an amount proportional to the 
interfering laser field y/K±. Vacuum fluctuations at frequencies close to the lower sideband frequency wo — uj m + Q only 
produce significant levels of radiation pressure noise by interacting with the lower sideband, and vice versa for the 
upper sideband. The coupling of the input field that is in-phase with the laser produces the term — kc\ in equation 
^ by an interaction of the fields that produces a displacement of the mirrors proportional to \JH and then a coupling 
of the laser field to the output through this displacement that is also proportional to ^fn. Generalizing to the case 
where laser fields at two different frequencies are resonating in the arms, the displacement of the mirrors is driven by 
radiation pressure noise from both frequency regions, and is thus proportional to ^JkZc\ + y/K + c^. This couples the 
field y/K— to the output port as k_ci + ^/k_k + C3 and the field to the output port as y/K-K+ci + K+C3, allowing 
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us to express the input-output relations as 



d 2 = Ad 2 + J2kZ-t — e i/3 , Ad 2 = (c 2 - k_ci - v ^r^ 4 Tc 3 )e 2i/3 
<fe = Ad 3 = c 3 e 2lP 



d d 



Adi + ^J2n+- — e , Adu = (04 — k + c3 — ^/K_K + ci)e 



^sql 



The relationship between the fields in the interferometer can be expressed by matrices as 
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e = pf + ra 



(7) 



Where 



b = —pa + rf 



(8) 



2/3 = 2arctan(fi/7) 



(9) 



is the phase shift acquired by a field at frequency £1 relative to the sideband frequencies upon reflecting from the arm 
cavities, 4> = (ujq ± u) m )l/c is the phase gained by the modulation sidebands while traveling one way in the signal 
cavity, and $ = fll/c is the additional phase gained by a (gravitational wave or quantum-noise induced) sideband at 
frequency uj ± u m + CI. These relationships allows us to solve for the output field b, which we express in compact 
form by 
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which is analogous to the notation used by Buonanno and Chen with the exception that our definition for the 
D-coefficients has the term y/2n± included in the coefficients rather than factored out which allows us to consider 
unequal intensities of the two sidebands. The compact form of our expression uses the definitions 
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D 3 = - (e 2l{ti+ ^p+ l) sin(0) (V^+^V^ sin(20) + (V^+^p 2 - e 2i ^+*) (2 cos(2<£) + sin(20) K _) p + lj ) (24) 
and 

£> 4 = - (e 2i ^+*V - l) cos(^) (e 2i ^+*VK- sin(20) + (e^+^V - e 2 ^+*> (2cos(2^) + sin(20) K _) p + l)) (25) 

These input output relations of equation 1101 form the basis for our analysis of the quantum-noise limited sensitivity. 
We note that by making the substitutions k_ — > k and k + — > these expressions reproduce equations 2.20-2.24 in Ref. 
for a signal-recycled interferometer with a conventional input laser spectrum. To investigate the quantum-noise 
limted sensitivity of an interferometer with this readout scheme we must first consider how these output fields are 
converted to a signal through heterodyne detection and demodulation. 



IV. DEMODULATION 



The output field of the interferometer is composed of two parts, a local oscillator and a signal (which includes the 
quantum noise) and can be written as 

E(t) = L(t) + S(t) (26) 

Where 

L(t) = Le Wot + h.c. (27) 

where h.c. stands for the hermitian conjugate, and L describes the amplitude and phase of the local oscillator. We 
follow the analysis of reference Q and express the signal as 

r+A jq i r 

S{t) = / = 6 1 (fi)e- l( "°- Wm+0)t +ib 3 (Sl)e- i( ' Ua - Um -W t + 6 3 (O)e- l(w0+ " m+o)t + ib^e"^"^^ + h.c 



+ (contributions at irrelevant frequency bands), (28) 

where A is the demodulation bandwidth. For simplicity we omit terms that are at frequencies that will not contribute 
to the final expression. 

The detected photocurrent is proportional to the square of the field giving 

»(*) ex E 2 out = L 2 (t) + 2L(t)S(t) + S 2 (t) (29) 

of which only the cross term 2L(t)S(t) has significant amplitude at the modulation frequency which is detected by 
demodulation and low pass filtering. This modulated component of the photocurrent is 

i mod {t) = 2L(t)S(t) = 2 (Le Wot + h.c.) x 

C+A JQ 1 

6i(fi)e" l(wo ^ m+n)t + i6 2 (fi)e" l( "°^ m+n)t + 6 3 (fi)e" l( " 0+ " m+o)t + i6 4 (f7)e~ l(a,0+Wm+f2)t + h.c. 
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(30) 



The demodulated output signal comes from mixing (i.e. multiplying) the modulated photocurrent with cos (uj m t + 4>d) 
and low pass filtering with a cut-off frequency of A giving 

cT 



idemod{t) = ■= I imod{t) COS (u! m t + <j>d)dt 
1 JT-t 
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where T — 1/A is the time constant associated with the low pass filter and we approximate the low pass filter as 
being a perfect step-function filter passing only frequencies within its bandwidth. Using the relations 
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we can rewrite the double-sided frequency integral in equation |^ as the single-sided integral 

idemod(t) » / ^e" 4nt [ftx^e 4 *- + iftafnje 4 *- + baffle-*' + ib^ffle^] + h.c. (34) 

Jo Zn V2 

which in the frequency domain is 

idemodity « -4 [6i (O)e^ + %(0)e^ +6 3 (n)e-*** +t6 4 (n)e-***] (35) 
V 2 

Setting C± = Co T 0d — tt/2 where Co = arg (£) we can express this as 

idemodity » ^ [i6i(r>)e iC - - 6 2 (fi) e lC - + z6 3 (r>)e- lC + - & 4 (fi)e-*+] (36) 
v2 

and with the definitions 

= ^3 sin C+ + &4 cos C+ (37) 
= b\ sin C_ + b 2 cos C- 

which represent the signal due to the output field in a band centered at loq ± uj m in the C± quadratures respectively, 
the demodulated signal is 

tde m od(^)~~^=(b ( -+b c+ ). (38) 

The values of bi in expression l38l are given by the input-output relations of ea uat ions 1101 and 1111 Since the output fields 
b^± each can contain vacuum fluctuations and signals from gravitational waves, the quantum-noise limited output 
is unlike that of equation 13 of reference [ 3]] for heterodyne readout of an interferometer with a conventional input 
spectrum for which the frequency regions around lu ± 2u m contribute noise (called non-stationary shot-noise) but 
not signal. Thus we have shown that heterodyne detection with a suppressed carrier input spectrum is not subject to 
non- stationary shot-noise. 

V. NOISE SPECTRAL DENSITY AND SENSITIVITY 

To compute the noise spectral density we express the magnitude of the noise by the equivalent signal level h n . From 
equation I1UI we have 



h n — hsql 



A&i sin C- + cos £_ + A63 sin C+ + A64 cos C+ 



t {y/2K-Di sin C- + y/2n-D 2 cos C- + \/2k + D 3 sin C+ + \f2n + D± cos C+) 
We find the noise spectral density by Equation. (22) of Ref. [jj 



(39) 



~ 2tt 6(n - n') s c h (f) = ~(m\h n (n) hi(n') + hi(n') h n (n)\m). (40) 



and assuming the input field a is in the vacuum state such that 



(0 a \a ia \,\0 a ) + (0 a \a\a f \0 a ) = ^2TrS(Q-n')5 ij , (41) 



the noise spectral density is 

S= - W + \A 2 ? + \A 3? + \A i? 

|t v /2k_( j Di sinC- +D 2 cos C-) + 2k+(D 3 sinC+ + D 2 cosC+)| 2 



The Ai terms in the numerator represent the noise contributions from each of the 4 input 'a' fields and are defined as 

Ai = Cu sin C- + C 2 i sin C- + C 3l sin C+ + C 4t sin C+ (43) 

We then define the sensitivity as the square root of the noise spectral density. For reference we reference this to the 
sensitivity associated with the standard quantum limit, h sq i 
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A. Balanced laser sidebands 



We now consider the special case of balanced laser sidebands such that k_ = k + . With this assumption C'n = 
C22 = C33 — C*44j C\2 — C34, C13 = C24 = C31 = C42, C14 = C32, C21 = C43, C23 = C41 Di = D 3 and 
Z?2 = -D4. With these simplifications, the noise spectral density is an even function of the demodulation phase 4>d an d 
is minimized for <pd — > 0. The noise spectral density is plotted in figure Eland compared to the noise spectral density 
with homodyne readout and with heterodyne readout using a conventional input spectrum. The noise spectral density 
for this heterodyne readout with a carrier-suppressed input spectrum is identical to that with a homodyne readout 
scheme and is quieter than that for heterodyne readout with a conventional input spectrum. The optimal value for 
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FIG. 3: The left hand figure shows the sensitivity for heterodyne readout with a conventional input spectrum (dashed line) 
and with the carrier-suppressed input spectrum (solid line). The sensitivity with the carrier-suppressed input spectrum exactly 
matches that achievable with homodyne detection (shown in the right hand figure). For the carrier-suppressed input spectrum 
I- = 1+ — I 3q i/2, while the conventional spectrum used for heterodyne and homodyne readout has Jo = Isql- For all spectra 
shown p = 0.9, r = y/l - p 2 , Q = 7/2, <j> = tt/2 - 0.47, Co = and <j) d = 0. 



Co the phase of the local oscillator is frequency dependent. The optimal sensitivity as a function of frequency and the 
local oscillator phase necessary for that sensitivity are shown in figure 0] 




FIG. 4 : Left : the optimal value of the local oscillator phase, Co as a function of frequency for an interfereomter with p = 0.9, 
r = yl — p 2 , /_=/+ = I sq i/2, fl = 7/2, (j> = 7r/2 — 0.47, and <f>d = 0. Right: the sensitivity at each frequency with the 
optimal local oscillation phase. 



B. Unbalanced laser sidebands 



In the general case of unbalanced laser sidebands, the lack of symmetry causes the optimal demodulation phase 
to differ from zero in a frequency dependent manner, however, unlike the phase of the local oscillator which is fixed, 
several different demodulation phases can be used to readout the signal in parallel channels, allowing for a noise 
spectral density that approximates that which would be achieved by having the optimal demodulation phase at each 
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frequency. The noise spectral density of the readout with unbalanced sidebands and frequency optimized demodulation 
angle is always between the best case of balanced sidebands and the worst case of totally unbalanced sidebands. The 
quantum noise ellipse @ is shown in figure El for the totally unbalanced case and the balanced case. 




0.1 I , , . , n/y "3 - 
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FIG. 5: Left: The sensitivity with unbalanced sidebands (7_ = 0.25I sq i, 1+ — 0.75I sq i) is plotted as the solid line. It is 
bounded by the sensitivity with balanced sidebands (I- — 0.5I sq i, 1+ = 0.5I 3q i) on the bottom (lower dashed line) and the 
sensitivity with completely imbalanced sidebands (7_ = 0I sq i, 1+ = Isql) on the top (upper dashed line). For all curves p = 0.9, 
r = y/l — p 2 , (f> — tt/2 — 0.47, and fo = 0. For the completely balanced and completely imbalanced sideband cases (f>d = 0, 
while for the unbalanced case the optimal demodulation phase, shown in the right panel is assumed at each frequency. 




FIG. 6: The total quantum noise ellipse (outer shaded ellipse) and its components from frequencies near the lower sideband 
(dark shaded inner ellipse) and the upper sideband (light shaded inner ellipse) . The ellipse representing a signal of magnitude 
h sq i is shown for comparison. On the left is the noise ellipse for unbalanced (I- = I sq i, 1+ = 0) sidebands on the right is the 
noise ellipse for balanced sidebands (J_ = J+ = I sq i/2). For both cases p = 0.9, r = \J\ — p 2 , SI = 7/2, <j> — tt/2 — 0.47, and 
4>d = 0. If homodyne detection were used the noise ellipse would be equivalent to the contribution from the lower frequency 
sideband in the unbalanced case (inner dark ellipse on the left) and the ellipse would be identical the total noise ellipse for 
balanced detection (outer ellipse on the right). 

It may seem that the ability to generate an output spectrum with a frequency optimized demodulation angle is a 
benefit of having unbalanced sidebands, however the optimal demodulation angle with balanced sidebands is frequency 
independent, so there is no disadvantage to balanced sidebands 



VI. CONCLUSION 



Future gravitational wave interferometric detectors will be limited by quantum-noise over much of the detection 
band. The details of how the readout scheme affects the quantum-noise limited sensitivity are an important consider- 
ation for the detector design . With the conventional input spectrum the quantum-noise with heterodyne readout is 
worse than with homodyne readout due to the presence of non-stationary shot noise - vacuum fluctuations at twice 
the modulation frequency which are down converted into the output signal by mixing with the modulation sidebands. 
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With the suppressed-carrier input spectrum presented here the quantum-noise limited sensitivity with heterodyne 
readout is equivalent to that with homodyne readout because the non-stationary shot noise is not present. Although 
both heterodyne readout schemes and homodyne readout schemes have other strengths and weaknesses, we show that 
the use of heterodyne detection does not necessarily lead to an increase in the quantum-limited noise floor. 
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